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WATER VAPOR LIDAR 

1 . Introduction and Summary 

This report summarizes a study of the feasibility of measuring 
atmospheric water vapor by means of a tunable lidar operated from the 
Space Shuttle. The purpose of such measurements would be to determine 
vertical profiles of humidity as functions of time, place and atmospheric 
activity. 

The specific method evaluated here is differential absorption , a 
"two-color" method in which the atmospheric path of interest is traversed 
by two laser pulses: one tuned exactly to the wavelength of a water 

vapor absorption line, the other tuned to a nearby wavelength in a gap 
in the water spectrum. Analysis of the two backscattered signals yields 
the number density [H2O3 as a function of range along the chosen atmo- 
spheric path. Of all the range-resolving lidar methods, differential 
absorption is the most sensitive for measuring 1H2O], EQ3] > [^02] 

and other molecular gases in the denser parts of the atmosphere. Thus 
the possibility of lidar monitoring of such atmospheric molecules from 
the Space Shuttle seems best evaluated from the standpoint of differential 
absorption. 

The question of feasibility can be phrased roughly as: Can we 

orbit a laser radar station of sufficient power and precision that useful 
atmospheric humidity profiles can be retrieved from the signals back- 
scattered to 200 kM altitude? The answer appears to be a qualified yes, 
depending on tradeoffs between laser pulse energy, repetition rate, 
range resolution, etc. The results are encouraging enough to warrant 
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a full-scale mission analysis specifying various orbits, percentage 
coverage of the atmosphere, and modes of operation which depend on cloud 
cover and the geographical regions of interest. 

The present report describes how one arrives at a qualified yes 
answer and touches only briefly on equations of mission design. The 
main topics dealt with here are: 

• I^hat are the needs of the basic and applied atmospheric 
sciences as regards humidity data, both locally and 
globally? (Section 2) 

• How are these needs met by conceivable lidar returns to 
the Space Shuttle, taking account of both atmospheric 
backscatter and H2O absorption? (Section 3) 

• What basis of confidence exists for sufficiently 
powerful, tunable and stable lasers that can be designed 
for the Space Shuttle within the bounds of present 
technology? (Section 4) 

It will be seen that, taken together, the answers to these questions 
imply that the overall concept is sound and deserves further study — 
not only from the viewpoint of mission design, but also demonstration of 
a breadboard system. The latter has Indeed been undertaken under a 
follow-on grant (NSG-1156) for joint work between NASA-Langley and the 
University of Maryland. Thus the content of Sections 3 and 4 relates 
partly to feasible, near-term studies from low altitude as well as to 
Space Shuttle performance per se . 

On the basis of the results reported in the following sections, we 
recommend that NASA undertake a mission definition study for Space Shuttle 
lidar measurements of atmospheric water vapor. For such a study the 
meteorological considerations of Section 2 will be paramount, and the 
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results of Sections 3 and 4 (and other, related studies by NASA-Langley) 
will be useful for scaling purposes. 
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2. A S^Tnmarv of Meteorological Requlremen ta for Water Vapor Data and 
Possible Space Shuttle Applications, 

2.1 Introductory Remarks 

The water content of air affects much of our daily life. Obviously, 
what we call weather - rain, snow, sleet and hail are direct consequences 
of the water content of air. In addition, the human and animal discom- 
fort indices and forest fire index require humidity observations as input 
variables. Industrial processes involving paper, paint, woodwork, textiles 
and precision Instruments are very sensitive to the surrounding humidity. 

In agriculture, harvesting and storage of grains require low atmospheric 
humidity. Excessive humidity has the dual effects on plant growth of 
changing evaporation from plant surfaces and of encouraging disease 
(e.g., Bourke, 1955 and Wallin, 1967). 

A considerable amount of information concerning the average seasonal, 
spatial and vertical distribution of the water vapor content of the tropo- 
sphere and lower stratosphere has been obtained from standard meteorologi- 
cal surface and upper air networks and from special stratospheric studies 
(e.g., Mastenbrook, 1966). Excellent average troposhperic water vapor 
profiles for different seasons and geographical locations are given by 
McClatchey et al. (1970) and the U. S. Standard Atmosphere (1966) and by 
Mastenbrook (1971) for the lower stratosphere. Climatological maps of 
surface and column water vapor content are given in Selected Climatic 
Maps of the United States (1965). Little climatological data on water 

vapor is available over the oceans. 

Unfortunately these cliTnstological data are not applicable to 
important atmospheric problems due to the fact that on an Individual 
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day, the water vapor density at a given troposhperic level may differ 
from the seasonal average by an order of magnitude, and the vertical 
variation is considerably different from day to day. This large variabil- 
ity is due to the t .ct that the water vapor distribution is a contributor 
to and a result of a given weather situation. That is, water vapor 
plays an important role in atmospheric processes at all timej and spatial 
scales of motion ranging from the molecular scale (i.e., evaporation) 
to the general circulation. Water vapor plays an important role in 
atmospheric energetics since it provides a source of energy through the 
release of latent heat, its presence affects the release of energy, it 
plays an important role in the radiation balance and, although a trace 
constitutent in the stratosphere, it plays a role in stratospheric photo- 
chemistry. 

Improved methods of measurement of water vapor alone will not lead 
to better theoretical models or weather forecasts. It is only through 
the simultaneous measurement of the thermodynamic and hydrodynamic pro- 
perties of the atmosphere combined with advanced knowledge and mathema- 
tical modeling of the processes governing the Interaction of the important 
variables that will result in improved forecasting. 

In this section we: 

1. Summarize the accuracy of water vapor measurement required by 
modelers and forecasters at a number of scales of motion. 

2. Review current and planned direct and indirect methods for 
operational use in obtaining atmospheric water vapor data. 

3. Identify important meteorological applications of water vapor 
data which might be obtained by a Space Shuttle Laboratory 


lidar system. 
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2.2 Water Vapor and Temperature Measurement Requirements 

Desired measurement accuracies for atmospheric variables for many 
scales of atmospheric motion have been specified for different experiments 
associated with the Global Atmospheric Research Programme (GARP) (see, 
Dods, 1970). Tlie measurement accuracy specification by different GARP 
committees does not encompass all possible scales of motion; however, 
they do cover the scales of motion which are possible to observe from 
the space shuttle. 

Fig. 1 is a diagram presented by Kuettner et. al. (1974) which 
illustrates thi. four major classifications of scales of atmospheric 
phenomena in the tropical atmosphere. Although the representation 
shown in Fig. 1 is only directly applicable to the tropics, the same 
scales of motion predominate in raid-and high latitudes. 

The largest scale, the A-scale, incorporates the so-called synop- 
tic and planetary scales. This scale is used for describing atmospheric 
waves of long lifetime, sometimes several weeks, and long wavelength. 

These waves include the so-called tropical easterly ';aves, Rossby- 
gravity waves of wavelength 5000 to 10000 km, and Kelvin waves. 

The next smaller scale, the B-scale is the scale on which tropical 
"cloud clusters" develop. It is thought that they form the link be- 
tween the short-lived smaller scale convective elements and the long- 
lived tropical waves and the Intertropical Convergence Zone (ITCZ). 

This scale of motion is generally not of great importance at higher 
latitudes (Kuettner et. al. , 1974). 

The C-scale refers to the scale of organized convection that form 
subsystems of the cloud clusters. In mid-latitudes, the term "mesoscale" 
is often used to denote this scale of phenomena. Two mid-latitude ex- 
amples of mesoscale phenomena are squall lines and groups of severe 


thunderstorms . 
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The smallest scale shown in Fig. 1, the D-scale, contains the 
individual convective elements; therefore, it is sometimes called the 
cumulus scale. This is also the scale of important variations of 
meteorological variables within an urban area. 

The accuracy and precision and time and spatial resolution of 
temperature and water vapor measurement necessary to satisfy the needs 
of forecasters and modelers of all meteorological processes at the 
scales of motion depicted in Fig. 1 have not been specified. However, 
the data requirements for some previous and planned GARP experiments 
have been cited in the literature. This section summarizes the data 
requirements for two such experiments for which the data requirements 
have been defined and which appear most stringent. These requirements 
should serve as a guide to the applications which might be made of data 
obtained from the Space Shuttle laboratory. 

The data requirements necessary for forecasters and modelers of 
large scale atmospheric phenomena have been thoroughly discussed by 
Dads (1970) in conjunction with planning for the First GARP Global 
Experiment (FGGE) to be conducted in 1979-1980. The primarv objective 
of FGGE is to provide initial and verifying conditions for experiments 
to enlarge predictability in the medium range (i.e., from a few davs 
to a few weeks), and in the extended range (i.e., from a few weeks to 
a few months). A comprehensive discussion of the FGGE objectives is 
given by the GARP Joint Organizing Committee (1974). 

The first major GARP experiment, the GARP Atlantic Tropical Experi- 
ment (GATE) , was conducted during the period June through September 
(1974) (see Kuettner et. al. (1974)). The GATE objective was to observe 
in sufficient detail the convective systems over the tropical Atlantic 
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Ocean, and to observe the larger scale atmospheric structure and motion 
systems in which the convective systems are embedded, so that the inter- 
action between the phenomena on these different scales is determined. 

An understanding of this interaction is to be used to develop para- 
meterization schemes for cumulus convection in large scale atmospheric 
models (Betts and Rodenhuis, 1974). The temperature and water vapor 
measurement requirements for the GATE Convection Subprogram at the 
cloud-cluster, meso-and cumulus scale were given by Betts and Rodenhuis 
(1974, pp. 59-63). 

Table 1 summarizes the large scale FGGE and smaller scale GATE 
temperature and water vapor measurement requirements. It should be 
noted that At refers to the time separation between simultaneous 
observations at each data point. Ax refers to the horizontal distance 
between successive data points and Az refers to the vertical separation 
of data required at each obsenration point. The accuracy requirements 
for temperature are in degrees Celcius and in percent relative humidity 
for water vapor. Appendix 1 discusses the relationship between relative 
humidity, water vapor mixing ratio, specific humidity and water vapor 
density. It should be noted that the stated accuracy requirements 
refer to relative accuracy. 

It is quite apparent from Table 1 that as the scale of motion de- 
creases, the required time, spatial and vertical resolution become 
smaller, and the accuracy requirements become more stringent. This is 
due to the fact that, in general, as the scale of phenomena become 
smaller, the shorter the lifetime, the more sensitive a model is to 
initial conditions (i.e., through finite differences), and the more 
rapidly observation errors propogate through a prediction scheme. 



TABLE 1. temperature and moisture observational requirements 


1 Scale 


RESoi^uxmN L 

AcoJi 

RACY 

At 

Ax 



^'ATER VAPOR 


Global 

12 HRS 

500 KM 

3 LEVELS IN THE ^ 
LAYERS lOOOiBOOv 

400-200, 125-SOmb 

Two LAYERS IN THE 
TROPOSPHERE 

0.5 - 1.0°c 

10 TO 20% 

Cloud-cluster 

100-1000 KM 

1 

6 HRS 

250 

Variable 
50 TO 200m 
from 

0 - 2 KM 

.5 - 3km 
from 

2 - 18 KM 

, 

0.2 TO 0.5 

5% 













TABLE 1* Continued 
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OLUTTON 

Accuracy I 
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AZ 

Temperature 

Water vapor 

Mesoscale 
10-100 KM 

i 

50 X 50 
100x100km 

1.5 - 

3 HRS 

10 KM 

Variable 

Highest re- 
solution near 
THE surface 
OF ORDER 

5Q-200M 

0.2°c 

2 - 5 % 

i 

1 

1 
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KM 

15 MIN 

1 KM 
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Highest resolu- 
tion NEAR the 
SURFACE OF ORDER 
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2.3 Current an d Planned Operational Attaospherlc Sounding Systems 
The largest amount of regularly obtained atmospheric soundings of 
temperature, water vapor, and wind velocity come from radiosonde ascents 
"from the worldwide network of upper air stations (see the Weather Station 
Index (1965) for a complete listing). These stations are irregularly 
spaced at roughly 400 to 500 km intervals in populated regions of the 
northern hemisphere but at much coarser intervals in South America, 

Africa and Asia. The oceanic regions are covered bv very widely spaced 
ship stations. 

Radiosondes are launched from each station at least twice dailv at 
0000 and 1200 GMT. Some stations also obtain data at 0600 and 1800 GMT. 
Temperature, water vapor and wind data are reported at a group of manda- 
tory pressure levels at approximately 1.5 to 2.0 km intervals between 
the surface and 30 km. Data are also reported at intermediate levels 
if significant deviations from standard soundings occur. The list of 
mandatory levels and criteria for significant levels is given in the 
Manual for the Radiosonde Code (1963). It should be noted that the 
designs and characteristics of radiosondes are not standard; therefore, 
some difficulties may arise when interpolating data across interna- 
tional boundaries. 

The radiosondes used by the United States military and National 
Weather Service (NWS) employ a thermistor and carbon strip hvgristor 
for temperature and humidity measurement. The accuracy of these sys- 
tems have undergone intensive review in recent years by Ostapoff et. al. 
(1970), Morrissey and Brousaides (1970), Teweles (1970), Lenhard (1970, 
1973), Betts et. al. (1974), and Brousaides (1975). The work of 
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Lenhard (1973) gives an upper bound of - 0.2 C for the 1)^^S radiosonde 
temperature error. This is well within the accuracy requirements for 
the previously discussed scales of motion. 

The accuracy of the NWS sonde humiditv device is less certain. The 
work of Ostapoff et. al. (1970) showed a diurnal relative humidity (RH) 
variation of 20% RH due to improper ducting and radiation errors. This 
led to a new design of the NWS radiosonde case (Friedman, 1972). 

-J- ^ 

For calibration purposes, the NWS requires a - 5% RH repeatability 
under all conditions except that below 30% RH, a - 7% RH error is al- 
lowed (Brousaides, 1975). Betts (1974) reports observed errors of - 5% RH 
at high RH values. However, Quiring (1973) and Brousaides (1975) show 
that errors much larger than - 9% RH occur when the humidity drops be- 
low 25%. As a result, the NWS no longer reports measured data less than 
20% RH. In addition, relative humidity data are terminated when the 
temperature is less than -40°C. 

To summarize, temperature and water vapor data from radiosonde- 
ascents are of sufficient quality for studies at scales of motion rang- 
ing from the meso-to the synoptic scale. However, mesoscale studies 
require a denser network than is available. Furthermore, the radiosonde 
gives poor humidity data at low humidities, and stratospheric water vapor 
data are not obtained. 

A large number of research and development and operational earth 
orbiting satellites, from which useful meteorological data may be ob- 
tained, have been launched or are in the planning stages. Table 2 
summarizes the reported capabilities of instrumentation on board the 
research and development satellites Nimbus 5, launched December 1972, 
and Nimbus F, launched June 1975, and the latest of an operational 
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series, ITOS-G, launched Noveniter 1974. Each of these satellites is 
placed in a sun-synchronous polar orbit. A comprehensive list of capa- 
bilities of previous and planned meteorological satellites is given by 
the Compendium of Meteorological Satellites and Instrumentation (1973). 

It appears that the present and planned remote sensing techniques 
are capable of providing water vapor data at the time and spatial resolu- 
tion required for large scale atmospheric studies. In addition, at least 
one technique will provide data on stratospheric water vapor. However, 
none of the remote sensing techniques give data with sufficient resolu- 
tion for application to mesoscale studies. 



TABLE 2. Reported Capabilities of Meteorological Satellite Systems 


1 

j Instrument 
1 

Range (km) 

Resolution (km) 
Ax Az 

Accuracy 1 

Temperature Humiditv 

1 Nimbus 5-ITPR 
Infrared 
Temperature 
Profile 
Radiometer 

0 to 25 

30 

3 to 5 

2.0°C 

■ 

0 to 6 

30 

6 


20% RH 

j Nimbus 5-SCR 
Selective 
Chopper 
Radiometer 

0 to 45 

30 

3 to 7 

2.0°C 


0 to 10 

30 

3 


20% RH 

Nimbus 5-NEMS 

Microwave 

Spectrometer 

0 to 18 

192 X 192 

5 

2.0°C 

1 

0 to 6 

192 X 192 

6 


1 

1 

Nimbus 5-THIR 

Temperature 

Humidity 

Infrared 

Radiometer 

5 to 10 

i 

L 

1 

5 

1 

— L 

i 

8 

i 

20% RH 


I 

M 

Ln 

I 
























TABLE 2 . Continued 


Instrument 

Range (fan) 

Resolution (km) 
Ax Az 

Accurac 

Temperature 

y 

Humidity 

ITOS-VTPR 

Vertical 

Temperature 

Profile 

Radiometer 

0 to 30 

50 X 50 

3 to 5 

1.0 to 2.0°C 


0 to 10 

50 X 50 

5 


20% RH 

Nimbus F-HIRS 
High Resolution 
Sounder 

0 to 45 

30 

2 to 7 

2.0°C 


0 to 12 

30 

4 


20% RH 

Nimbus F-PM2 
Pressure 
todul ated 
Radiometer 

0 to 65 

71 to 360 

3 to 5 

5.0°C absolute 
2.0°C relative 

- 






Nimbus F-SCAMS 
Spectral Scanning 
Microwave 
Spectrometer 










: 

Nimbus F-LRIR 
Limb Radiance 

1 

Infrared 

Radiometer 

lA to 75 

400 

2 to 3 

3^C 

■ 

14 to 50 

400 

2 to 3 


mixing 
20% ratio 
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2.4 Applications of Water Vapor Profiles Obtained from the Space 
Shuttle Lldar 

The meteorological applications which might be made of water vapor 
profiles obtained with the Space Shuttle lidar system depend priraarilv 
upon : 

1. the accuracy and pfe;crsiSTr of the Ildar retrieval technique. 

2. the horizontal and vertical resolution possible from Space Shuttle 
altitudes. 

3. the inclination and duration of the Space Shuttle orbit. 

In order to have application to all the scales of motion shown in 
Fig. 1, the retrieval technique must meet the requirements of the cumulus 
scale. If one uses the rough estimate of 10 km sec ^ for the ground speed 
of the shuttle, one would require a 100 hertz firing rate of the lidar 
system in conjunction with a mirror system in order to obtain one snap- 
shot of Che atmosphere at the horizontal resolution required at the 
cumulus scale. The power requirements for that firing rate appear to 
be beyond the present or near future state of the art. However, it 
appears that 10 hertz firing rate may be attainable for short durations 
by the time of Space Shuttle flights. Therefore, with a proper mirror 
system, the smallest scale the lidar might be applied to is the raesoscale. 
In addition, such a firing rate could provide vertical profiles at the 
horizontal resolution of the cumulus scale. Assuming the accuracy and 
vertical resolution requirements of Table 1 can be met, the Space Shut- 
tle lidar can provide water vapor data with some applicability to all 
the scales of motion shown in Fig. 1. 

The application of water vapor data to studies of large scale atmos- 
pheric processes depends primarily upon the inclination and duration of 
the orbit. Obviously, a low inclination orbit around the same latitude 
circle could provide, at best, a time series of data from data 


sparce 
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oceaaic and unpopulated regions. On the other hand, a polar orbit would 
provide global coverage, although the data would be asynoptic in character 
since they would not be gathered at one time. However, four dimensional 
data assimilation techniques (e.g., Rutherford, 1972; Bengtsson and 
Gustavsson, 1972) would make proper use of these data. In either case, 
however, it appears that current or planned passive satellite retrieval 
techniques can provide water vapor data at the resolution required for 
large scale atmospheric studies. 

It appears that the greatest application of Space Shuttle Ildar 
data would be to studies and modeling of phenomenon at the cloud-cluster 
and meso— scales. A number of such studies have been performed in recent 
years and several are planned for the future. These studies include 
GATE-1974 (cloud-cluster and mesoscale) , the National Hail Research 
Experiment (NHRE) - 1971-76 (mesoscale) - Foote and Prankhauser (1973), 
the Air Mass Transformation Experiment (AMTEX) - 1974-75 (synoptic 
and mesoscale) - Lenschow and Agee (1974), the Severe Environmental 
Storms and Mesoscale Experiment (SESAME) 1978-80, see Lilly (1975), 
and the Polar Experiment (POLEX) - 1975-76, (synoptic-mesoscale) , 

Weller and Bierly (1973). 

These experiments require the simultaneous measurement of many 
atmospheric variables, including temperature and water vapor at meso- 
scale resolution over a limited area larger than the cloud-cluster scale. 
Fig. 2 shows the horizontal surface distribution of several meteorological 
parameters determined by Foote and Frankhauser (1973) in the vicinity 
of a mid-latitude mesoscale convective feature. This figure points out 
the complexity of the horizontal distribution of several^ important meteoro- 
logical parameters and illustrates the necessity for measurements at 
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small horizontal spacing. It should be noted that a detailed three dimen- 
sional mapping of one parameter, such as the water vapor mixing ratio 
(Fig. 2(e)), is not sufficient to predict the location of the greatest 
water vapor convergence (Fig. 2(c)) or the time variation of a mesoscale 
convective feature. It is apparent that the usefulness of a Space Shuttle 
water vapor lidar to studies of mesoscale convection would be providing 
supplemental data to an existing field experiment. 

An important meteorological feature, in which water vapor plays 
an important role, is the so-called "dryline" front. The dryline has been 
discussed in detail by Rhea (1966) and Schaefer (1974a, 1974b, 1974c). 
Briefly, the dryline is a narrow zone across which a sharp moisture 
gradient occurs (roughly a factor of two change in mixing ratio over 10 Ion). 
It often exists and moves under synoptically quiescent conditions; how- 
ever, it is often associated with thunderstorm development over the southern 
plains. The dryline is oriented north-south, approximately parallels 
the terrain contours and is often found over west Texas and Oklahoma 
during the spring and early summer. The dryline surface has a character- 
istic steep profile, being roughly vertical to approximately 4000 ft. 
above its surface intercept. It then tilts rapidly eastward to become 
nearly a horizontal surface. Drylines typically move roughly 500 km 
eastward during the day and westward at night. Although narrow, the dry- 
line may be over 500 km in length- According to Schaefer (1974c), the 
dryline is also a significant feature during the pre-monsoon months over 
India and Central West Africa. 

Schaefer (1974a) has developed a numerical model of dryline move- 
ment. One large uncertainty in forecasting its position is its initial 
position. Currently, its position is obtainable only after it has passed 
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one of the widely spaced surface stations. According to Schaefer (1974b), 
the temperature does not vary across the dryline. Therefore, the factor 
of two change in the mixing ratio across the dryline Is due only to a 
variation in the water vapor density. If the Space Shuttle water vapor 
lidar could obtain the surface to 4000 ft column content with less than 
ten percent random error, observations at one to two kilometer horizontal 

resolution could be used to position and track the dryline front on each 
orbital pass. 

A 10 hertz Ildar sampling rate would allow a cross section of verti- 
cal water va^or profiles with roughly one kilometer horizontal resolution. 
Such data could be used to sudy the variation of water vapor in urban 
areas and urban effects on the larger scale water vapor field. Experiments 
of this kind would be useful in delineating areas conducive to fog and 
smog formation and downstream convective development. 

A final application, which has not been studied in detail, is the 
determination of the stratospheric water vapor content. The only current 
passive instrument suited for this purpose is the Nimbus F Limb Radiance 
Infrared Radiometer (LRIR) which claims ar. accuracy of +20% mixing ratio at 
2 to 3 km vertical resolution for the 15 to 50 km range. In order to be 
competitive, the Space Shuttle water vapor Ildar would require equal or 
better accuracy. 

Perhaps the bast method for obtaining stratospheric water vapor 
iniormatlon from a Space Shuttle lidar would be the measurement of the 
backscatter from the tops of optically thick cirrus. If these observations 
could be compared to similar observations in clear areas, one might obtain 
estimates of the amount of water vapor Injected into the stratosphere by 
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thunderstorms. However, it is likely that the clear column lidar obser- 
vations would be dominated by tropospheric water vapor. Therefore, it 
would be extremely difficult to extract the injected amount with the lidar 
unless an independent observation were obtained. The applicability of 
a Space Shuttle water vapor lidar to this problem requires further simula- 
tion using estimates of the stratospheric aerosol distribution and the 
optical properties of cirrus clouds. 

To summarize. Space Shuttle lidar determined water vapor profiles 
might be applied to; 

1. augmenting synoptic observations in data sparse areas. 

2. three-dimensional mapping of mesoscale water vapor fields. 

3. detection and tracking of "dryline" fronts. 

4. determination of urban influences on the larger scale water 
vapor field. 
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N76“ 18482 

3. Lidar Simulation 

Tills section estimates the feasibility of measuring atmospheric 
water vapor via orbital lidar. The calculation starts with laser radar 
equations representing backscatter with and without molecular line ab- 
sorption; the magnitudes of off-line backscatter are demonstrated (3.1). 
Extensive prior data on water line strengths are summarized in 3.2 to 
indicate the available sensitivity to water vapor concentration. The 

O 

wavelength range of greatest interest is 7000 A - l.Oy. This range is 
accessible to various tunable dye laser techniques, such as ruby laser 
pumping of polymethine dye solutions in a well— designed, narrow band 
cavity (see Section 4). 

Several lidar situations are then considered (3.3) starting with 
uniform and perturbed atmospheres at 0, 3, 10 and 20 kM (stratosphere) 
altitudes. These simulations are indicative of results to be obtained in 
"ground truth" measurements (ground-based and airborne) . An approximate 
treatment of polar observations is also given. Finally, vertical atmosi* 
pheric soundings from orbit and from ground stations are calculated. 

Errors are discussed (3.4) as regards their propagation through the lidar 
equation to render the measured water vapor concentration imprecise; 
conclusions are given as to required laser energy and feasible altitude 
resolution. 


3 . 1 Prelimina r ie s 

We adopt a monochromatic lidar equation for the ratio 
between the received and transmitted numbers of photons in the laser pulses: 


R 



3/r) 


r 

A 1 -2 I 


( 1 ) 
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where K is a total optical system efficiency, t is the "range gate 
or open time for accumulating is the volume coefficient of 

backscatter and is in general a function of the lidar range r, and 
A/r^ is the element of solid angle subtended by the receiving telescope 
of area A. Since cx/2 is the range resolution (kM) , and the dimensions 
of 3 and A/r^ are respectively (km- star ad) and (sterad) , the 

expression in brackets is dimensionless. 

The last term represents all extinction processes during two-way 
passage of the beam along the line of sight. The total extinction per 
unit length can be decomposed into absorption and scattering by molecules 

and particles: 


8=6 +6 + 

*^T a,m a,p s,m s,p 


( 2 ) 


-1 , 


= S + III (all in units kM ). 
a,m 


1-Iolecular absorption is singled out for attention because it is usually 
much more dependent on wavelength than the other three terms (III) . 
Explicit use of this property is mads below in seating up the simplest 
form of the differential absorption method. 

The lidar equation (1) assumes no refractive or multiple scattering 
effects in the atmosphere; also the detector must be so apertured and 
gated that at each range it receives all the light, and only the light, 
backscattered from that conical slice through the atmosphere which is 
illuminated by the slightly divergent laser beam. 


Scattering processes 


therefore enter twice into equation (1);— they 
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uitimately determine the 180“ backscatter coefficient 8 and give rise 

7 T 

to part of the extinction via 8 and 8 . To simplify calculations, 

s,m 5,p 

we will let 


I IT hi 

i ir 


(3) 


where e == 1 for molecules and small aerosols and g = 0.2 for large aero- 
sols (Mie scattering). The latter reflects the fact that roughly 80" 
of the intensity of Mie scattering lies in the forward direction, leaving 
20% for distribution over all angles. We assume that the rich fine struc- 
ture of Mie scattering is usually washed out by distributions of particle 
size and composition found in nature. 

Since these scattering coefficients and the particulate absorption 

8 _ are only weakly varying functions of the wavelength, we can set them 

a,p 

equal in two lidar equations for closely spaced wavelengths (X^, X) 
over which molecular absorption may vary drastically. X^ might be at 
the center of an H^O line, for example, while X might be 3 or 10 half- 

O 

widths away. Since | X^-X | can therefore be as small as 1-2 A, e.g. for 

H 2 G and O 2 lines in the near infrared, it is reasonable to focus attention 

on 8 (X) and allow constancy of the other extinction coefficients over 

a ,m 

The ratio of the lidar equations for the two wavelengths is then 


Q(r) “ R (r)/R (r) = exp 2 f (8 (r’) - 8 (r’))dr', 

A A I &1TU1A 

ri J n ’’ 


(A) 


and (1/Q) (d/dr) (Q) yields 



-26- 


^ [In Q(r)] = 2 [6^ (r) - 6^(r)] . 

o 


( 5 ) 


In general the molecular extinction will be 6 = n{r) o(r), where both 
the number density and the absorption cross section may vary along the 
line of sight. For sufficiently constant temperature and total pressure, 
changes in molecular line shape may be ignored, B(r) = o (const.) n(r), and 


^ [In Q(r)] 

" 2[a(X^)-o(A)] • 


( 6 ) 


This is the central equation of the differential absorption method whereby 

one obtains range-resolved n(r) from the observed time (i.e., distance) 

behavior of the ratio Q between off-line and on-line lidar returns 

(Measures and Pilon (1972), Ahmed (1973), Ryer and Garbuny (1973). 

For best results the normalized off-line return should be 

as large as possible, so that modest amounts of absorption in the on-line 

return are discernible with good statistical reliability. Essentially 

one is relying on atmospheric backscatter as a "distributed mirror" along 

the line of sight in order to get sufficient photons returned from any 

given layer without overly attenuating the beam for layers further away. 

For a given range, the atmospheric scattering conditions 6 giving 

s 

the optimum off-line return from that range can be calculated from equa- 
tions (1) - (3). Since 3 and 6 vary with wavelength over inter- 

s y m ^ 9 

a 

vals of hundreds and thousands of Angstroms, the optimal lidar wavelength 
for a given range will also depend on atmospheric conditions. Representa- 
tive values of 3 are given in Table 3 which is taken largely from 
s 

McClatchey et al (1972). 

i 

*Note added in proof: References to R. M. Schotland’s pioneering and 

continuing work (1966, 1974) were inadvertently omitted. Tliey are given 
on p . 86 . 
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The principle of optimum B for a given range can be illustrated 
by assuming a uniform atmosphere in which one type of scattering is domin- 
ant. Then equation (1) can be written 


R = C 



whence “ 0 yields l-2rB = 0, or 


(7) 



e.g. to observe optimum lidar returns from 5 kM, one should choose condi- 
tions such that 6 = 0.1 kn“^; then also the return from 4 kM will in- 
s 

crease wa'.th increasing haziness while that from 6 kM range will decrease. 

These considerations are important in the design of ground-based 
experiments where the weather conditions are variable and where one would 
try to validate system errors prior to flying a system in an aircraft 
or satellite. Remsberg et al (1975) have recently reported estimates of 
the limits of utility of orbital , differential absorption lidars based on 
similar considerations of the off-line lidar return. 

To illustrate the lidar returns for various situations, we plot 
In Figure 3 the off-line returns for four wavelengths and two uniform 
atmosphere situations according to McClatchey's (1972) summary of data 
on typical conditions. The lidar parameters are given in the figure. 

The curves were calculated according to equations (1) - (3). As to the 
absolute numbers of detected photons, it is noteworthy that the assumed 
15 meter range resolution \ 30 nsec range gate) could well be Increased 



Table 3 , Coefficients (kM 

by Scattering in the Near 

for Optical Extinction 
Infrared (A=8600 A). 

** 

Model 

6s, molec. 

Bs, particle 

6tiot=Bs,m +6s,p 

Hazy, midiatltude summer (sea level) 

0 . 00193 

0.440 

0.442 

Clear, midlatitade summer (sea level) 

0.00193 

0.0903 

0.0922 

^ ^ ^ 

Antarctic summer (3kM ground level) 

0.00145 

0.00480 

0.00625 

Clear, midlatitude summer (alt. 3kM) 

0.00143 

0.00817 

0.00960 

Clear, tropical (alt. lOkM) 

0.000677 

0.00179 

0.00247 

Tropical (alt. 20kM) 

0.000154 

0.000860 

0.00101 

* « 

^ = 6; 8 = on 




** Most data interpolated from tables 
R, A. McClatchey, R. W. Fenn, J. E. 
Report no. 72-0497, (August, 1972). 

glnen in "©ptical 
A. Selbv, F. E. 

Properties of 
Volz, and J. S. 

the Atmosphere", 
Caring, AFCRL 


*** Molecular scattering estimated as comparable to Subarctic summer (**) : Particulate 
scattering based on N = 150/cm (W. Zoller, private communication) and a =3.2 x 10“^® cm^ 

estimated from AFCRL** values for B vs N. P 

s p — — 



# photons detected /1 50 M range cell 
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LIDAR RETURNS (near IR) 

due to bacKscatter and extinction 

by aerosols at sea level 

z 

300 mJ pulse, !M rcvr. 

I^sec range cell 


Clear atmosphere 
( 23 kM visibility) 


-11,000 A 
,9400 
,8200 
,7300 A 


•11,000 A 
■S400 A 
-8200 A 
■7300 A 


range (kM) 


Figure 3 


3>o 3>o 
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for larger scale observations at the longer ranges so as to minimize 
photon noise; also much smaller range resolutions, say IM corresponding 
to the 6 nsec duration of a N^-laser-pumped dye, could be used out to 
ranges of a few kM without incurring too great a photon noise. More 
examples of off-line Ildar returns will be shown below as part of the 
discussion of the detectability of water vapor. 

A further conclusion we make from Figure 3 is that the variation 
between models (for this X-range) is much greater than the dependence 
on wavelength within a given model. Henceforth for any given aerosol 
m odel we will employ a s tandard calculating wavelength of 8600 A to 

represent _ scattering in an a v erage way for all H.,0 lines b etween 7200 A 
and 1.0 u . 

This facilitates a "parametric" view in which one assumes with rea- 
sonable justification that more-or-less constant tunable laser performance, 
detector efficiency, and atmospheric backscatter (for a given aerosol) 
can be obtained over this near infrared region — and that the principal 
difference from case to case will reside in the H^O absorption cross 
sections ("line strengths") available in any given band, and in the 
abundance of water vapor and aerosols in the atmosphere. 
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3,2 Water Spectrum; Wavelengths, Strengths and Widths . 

This section arrives at useful representations of the strengths of 
water lines available in the visible and near infrared. Cross sections 
are indicated for the subset of lines we will call "temperature insensi- 
tive" — namely those lines whose lower state populations are stationary 
with temperature for temperatures in the vicinity of 300“K. This will 
guarantee that the cross sections are representative of the values to 
be employed in atmospheric observations, since temperature— insensitivity 
will be an important property of any procedure for determining molecular 
abundance . 

The extent of the H^O absorption spectrum into the visible is 
indicated in Fig. 4, which is a low resolution scan of atmospheric trans- 
mission (sea level) over a 300 M path. Tunable dye laser techniques 
enable wavelengths up to 1.1 p to be used for absorption measurements, and 
the parametric oscillator technique works well beyond this wavelength. 

What concerns us then is the strength of lines in the various bands such 
as at 7200 A, 8300 A and 9400 A. We will see that many of these lines are 
sufficiently strong for [H 2 O] assay from orbit and from the ground, 
while the 1.4, 1.9, 2.7y, ... bands are too strong unless one is concerned 

O 

with stratospheric H 2 O. The lines near 5900 A are interesting because of 
their coincidence with the region of peak performance by dye lasers 
(Rhodamine 6-G) , but fall on the weak side as regards general applicabi- 
lity to atmospheric problems. 

To illustrate the orders of magnitude involved, recent measurements 
by Meredith et al (1973) at Science Applications Inc. yield peak cross 




Atmospheric 
alcen from ! 

Technology (Ed. W. L. Wolfe) Office of Naval Research and 
U. S. Govt. Printing Office (1965); original source: H. W, 

Yates and J. H. Taylor, "Infrared Transmission of the 
Atmosphere", NRL Report 5453, ASTIA AD 240188 (I960).] 
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sections of order 4 x 10 cm for strong, temperature insensitive lines 

O 

of H 2 O near 9400 A. Thus the 1/e attenuation length (optical depth 

unity) for a laser tuned to such line centers is of order 100 M under 

tjFpical ambient conditions of humidity (>2.7 x 10^^/cm^, or p(H20)=10~^atmos. ) 

For atmospheric ranging over total path lengths of 10 kM, such lines will 

clearly be of greatest interest in the polar regions where the absolute 

humidity is of order 1/100 of the above. For work in the temperate and 

tropical zones, one needs H 2 O cross sections in the range 2 x lO'"^^ - 
-23 2 

2 X 10 cm . Indeed these are to be found In the shorter wavelength group- 
ings shown in Figure 4. 

For quantitative discussion we will adopt the standard, serai-empirical 
relationship between absorption per particle in the line of sight and the 
strength and width of a Lorentzian line; 


o(v) 


1 Sy 

{{v-v ) +Y } 
o 


(9) 


o is an absorption cross section (cm^) ; S = |o(v)dv is the line strength 
(cm cm ); Y is the line half-halfwldth (cm ^), For absorption in a 
uniform gas of length £, the extinction of intensity I(v) is given by 


I(v)/I^(v) = = ^-a(v)n£ ^ 


( 10 ) 


where n is the total number density of the species in question (not the 
number in a given state). In this representation, -he line strength S 
clearly involves the Boltzmann factor for the lower level of the line as 
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vell as ^he partition function for the species. S is sometimes specified 
in units (cm in applications where one knows the mass density p 

and the total mass per unit area (p&) in the line of sight. 

Xn a previous report Wilkerson, Ercoli and Tomkins (1974) gave 
detectabilities" for H^O in terms of total gra/cm , by requiring an opti- 
cal depth 6(v) = 0,03 and dividing by typical strong line strengths in 
various wavelength regions. The present summary of cross sections is more 
quantitative since it is based largely on new data which are extensive 
enough that the question of terapei ature-insensitivity can be addressed ex- 
plicitly. We are greatly indebted to Professor W. S. Benedict (University 
of Maryland) and Dr. J. W. Brault (Kitt Peak Obs.) for the opportunity of 
using these data. Other published sources are also used for selected longer 
wavelength lines. 

All the lines comprising the H^O bands in the upper part of Figure 4 
involve transitions to high vibrational states (v^^, v^] of the electronic 

ground state of H^O; e.g. {401} near 5900 A, {301} near 7200 A, and {211} 

O 

near 8200 A. What is of greatest interest from the standpoint of appli- 
cations is the subset of lines within each band which originate on rota- 
tional sublevels whose populations n^^ are temperature-stationary; i.e., 

3 ^ — 0 . ( 11 ) 

Restricting one's attention to these lines, and requiring that they be 
easily resolval :.s from all other lines in the H^O spectrum, guarantees that 
one's lidar feasibility estimates will be based entirely on useful absorp- 


tion lines. 
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Table 4 shows rotational levels for {o, o, o} and other vibrational 

states of H 2 O, taken from Swensson et al (1970) . The levels singled out 

are those for which the energy is about 300 cm ^ above the rotational ground 

state. In descending order, these may also be written in (J", K ", K ") 

3 c 

notation as 3 ^q, 3^^^, ^13’ ^23’ ^22’ ^05’ ^15’ total 

angular momentum and K. , K are the components relative to the axes of least 

3 C 

and greatest moments of inertia of the H^O molecule. 

The indicated levels have nearly stationary populations, relative to 

3/2 

T, for T = 300°K. This arises from the dominant T dependence in the H^O 
rotational partitiion function Zj^(T) and the level population in each lower 
rotational level 




K ” 
c 


ne 


-E"/e 


z^(0) 


( 12 ) 


where n is the total number density of H^O molecules, E" is the rotational 
energy (cm ^), and we have introduced the corresponding temperature scale 
(cm ^ units) 0 = kT/hc. Equations (11) and (12) yield the level-energy of 
maximum population 


E" 

max 



i.e., levels having E" 'Xi 310 cm ^ will have stationary populations for 

T = 300°K. Thus the seven levels indicated in Table 4 are the most suitable 
for temperature-insensitive determinations of H^O abundance by means of 


F. Mills (NASA-Langley) has pointed out that lower levels (200-225cm 
in J = 3,4) may be better, if the laser line can be made very narrow compared 
to absojrption lines, owing to the 'X 1 TV 2 dependence of the line width y which 
offsets part of the T ^/2 dependence in the molecular rotational partition 
function. 
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Tabla 4 [taken from The Solar Spectrum from A7wS*8 to X120I6 by J, W. 

Swenaaon, W. S. Benedict, L.Delboullle and G. Roland, Instltut d'Aatrophyslque 
de I'Univeralte de Liege (1970).] Levela near 300 cm~* are Indicated, 
having 3rij^/3T«o for T”300®K. 

HsO Energy Levels (cm-i) 


Noit I L»vela fully confimiBd by combination dilTorenoes ore given to 0.01 cm"*. TTnoonfirmod lovaU of high pm- 
bobility orogtven toO.l cm'*, with the mom doubtful indicatol by t. LowuIb m difforant vibrational ataua whoso pruportii H 
ait) roixed by perturbation to on aiceptionnl degree are denoted by 

Section 1 


J 

T 

000 

130 

031 

210 

111 

012 

0 

0 

0.00 

8273.95 

8373.82 

8781.57 

8807.00 

0000.13 

1 

-1 

23.79 

8297,34 

8307.46 

8784.00 

8830.21 

9023.47 


0 

37.13 

8323.31 

8421.18 

8799.65 

8844.47 

0037.21 


1 

42.37 

8329.38 

8427.34 

8805.14 

8850.00 

0042.31 

2 

-2 

70.08 

8343.14 

8443.63 

8829.48 

8876.30 

0008.70 


-1 

79.48 

8364.05 

8402.44 

8840.27 

8886.10 

9078.29 


0 

96.17 

8382.21 

8480.80 

8856.75 

8903,46 * 

9095.10 


I 

134.91 

8468.03 

8640.96 

8899.83 • 

8944.67 

9135.34 


2 

136.17 

8467.78 

8650.90 

8902.60 

8946.08 

9137.25 

3 

-3 

136.77 

8409.66 

8610.34 

8894.00 

8039.93 

9133.60 


-8 

142.28 

8424.60 

8623.08 

8900.66 

8045.64 

91,30.04 


-1 

173.36 

8480.77 

8660.12 

8933.19 

8979.06 

0172.49 


0 

208.30 

8627.32 

8620.92 

8909.72 

9013.98 

0205.00 


1 

212.16 

8532.08 

8625.98 

8976.67 

9020.17 

9212 46 


2 

286,23 


8660.62 

8746.62 

9067.78 

9096.30 

0287.18 


3 

286.43 


8660.71 

8746.64 

9067.88 

0098.61 

0287.44 


-4 

222.08 

8494.67 

8606.54 

8976.29 

9022.30 

9216 13 


-3 

224.83 

8504.35 

8603.50 

8079.63 

9026.15 

9218 34 



276.52 


8663.87 

8884.07 

9032.09 

0070.71 

9273.26 


-1 

300.36 ] 


8620.23 

8714.59 

9061.13 

9106,42 

9298.21 


0 

315.70 


8633.10 

8728.46 

9077.32 

9121.49 

0315.17 


1 

382.52 

8760.62 

8843.23 

9162.60 

9103.15 

0383.00 


2 

383.85 

8767.11 

BS44.03 

9163.62 

0194.63 

9384.07 


3 

488.10 

8929.44 

9006.65 

9207.96 

9305.40 

0400.94 


4 

488.13 

8929.46 

9006.68 

0207.98 

9306.48 

0400.07 

0 

\-6 

326.36 


8696.88 

8697.99 

9075.00 

9121,72 

0310.10* 



326.64 


8602.77 

8702.46 

9077.16 

0123.01 

0317 03 



399.46 

8690.08 

8701.02 

9164.16 

0200.05 

0304 00 


-8 

416.23 

8736.05 

8830.14 

9173.37 

0218.02 

0411.04 


-i 

440.60 

8761.65 

8868.60 

0204.67 

0240.23 • 

0444.08 


0 

604.00 

8876.35 

8064.04 

0280.00 

9310.08 • 

0602.71 


1 

608.81 

8878.58 

8007.10 

0276.03 

03IG..10 

0603.21 


s 

610.12 

0050.02 

9127.06 

03SG.7 

0424.37 

9611.61 


3 

610.36 

0061.02 

0128.00 

9380.06 

0424.64 

0011.70 


4 

742.10 

0267,3.5 « 

0324.00 

0628.51 

9.563.07 

0746.18 


S 

1 742.10 

0267.18 

0324.13 

0528.61 

0603.07 

0740.18 


■ of the 
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optical absorption. 

Lines originating on these levels are to be found in all the vibra- 
tion-rotation bands shown in Figure 4. We cuJled available data on line 
strengths S and halfwidths y within these bands, to establish reasonable 
limits on the absorption cross sections for temperature- insensitive lines 
in each wavelength region. The so-derived limits on £ as a parameter 
were then employed in lidar simulations. The remainder of this section 
deals with the line strength results. 

9400 A ; {201}, {300} 

The results of Meredith. et al (1973) may be summarized by an absorp- 
tion cross section = S/vy in the range 1.5 - 6,5 x 10~^^ cm^ for several 
prominent lines which are temperature-insensitive or nearly so. 

O 

8200 A ; principally {211} 

Farmer's (1971) study of the stronger lines in this region gave strengths 
in the range (0.6 - 2.5) x 10 cm for seven temperature- insen- 

sitive lines. Assuming a common FWHM = 2y = 0.16 c.a~^ for an STP atmos- 
phere, one obtains = S/ry in the range (2.5 - 10) x lO"^^ cm^ for typical 
strong lines that would be employed for lidar measurements. Much smaller 
cross sections are also to be found in any such region, as can readily be 
judged from the large spread of H2O absorption intensities seen in the 
Fraunhofer spectrum. 
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6810 - 7470 A ; {301, 103. 221, 202, 400, 023, 004, 122, 141, 320} 

The range of line strengths for 90/100 temperature- insensitive lines 
terminating in these upper vibrational states can be adequately repre- 
sented by S 3x10 - 1.7x10 cm ' » according to the data of 

Brault and Benedict. Their listing covers roughly 2250 lines in this region, 
clustered mainly in XA 7000-7300. For Yg.pp ''' 0.09 cm the corresponding 
range of peak absorption cross section is: 


fa = 1 X 10 - 6 X 10 

[^for AA 6810-7470 (H 2 O) 


(13) 


Thus the strongest cross sections compare with ^hose given above for the 

O 

8200 A region. The magnitudes of the cross sections are distributed 
fairly avanly over this range, so that the factor-of-500 is reasonable to 
employ as a parametric variation of a^ in lidar calculations. 


6270-6660 A : (231, 212, 311, 113, 410} 


Data on about 1000 H 2 O lines in this region were provided by Brault 
and Benedict, Of these, 55 are useful lines originating on levels whose 


The water line at y 6943.8060 is well known for its proximity to the 
ruby laser emission. The two can be made to coincide by temperature- 
tuning the ruby, so that this line is a natural candidate for lidar measure- 
ments of H 2 O. The line strength is given as 0.1773 cm“V(gm/cm^) by Brault 
and Benedict. Moreover the line arises from absorption out of the 4^:, level 
of {0,0,0}, designated in Table II as J=4, t=- 3, E“224.83 cm“^. In Mills' 
interpretation alluded to above, this line would show a temperature-insensi- 
tive central absorptivity for laser radiation having a llnewidth « 0.1 cm~^. 

This leads us to provide an estimated effective cross section S/ny here 
for the sake of completeness: A 6943.8060 A ; 5 = 5. 31xl0“^’'cm"^/(molec./cm2) ; 

YgTp “ 0.098 cm~^i o^ »= 1, 74x10"^ ^cm^ . One obtains an optical depth 6 “ "^ 0 *^ “ ^ 
for i “ 2.1 IcM and n ■= 2.7x10^^ cm”^ (typical 1% absolute humidity). 
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populations are temperature- Insensitive; moreover, 50/55 have strengths in 

the range 1.5 x 10 - 1.8 x 10 cm For Yqt,p " 0.08 cm 

cm olr 

one obtains 


fa = 6 X 10 - 7.2 X 10 cm' 

\ ° 

[_for AX 6270 - 6660; 


(U) 


i.e., the stronger lines here have cross sections which are about an 

O 

order of magnitude lower than for the strong lines around 7200 and 8200 A. 
Summary (3.2) 

For the near Infrared below 1 micron wavelength, the effective ab- 

”*26 2 

sorption cross sections for H 2 O lie mainly in the range 6x10 cm to 
-22 2 

6x10 cm , considering only the temperature-insensitive levels. For many 

of the calculations below, the parametric range for was chosen to be 
2A *^21 2 

10~ - 10 '*‘cm , since weaker lines are not useful and stronger H 2 O lines 

are not to be found in this spectral region. 

3.3 Results of Lidar Calculations 

This subsection presents the results of li' r calculations for a 
variety of atmospheric situations : 

3.3.1 Horizontal ranging at various altitudes in atmospheres of 
uniform composition (including stratospheric and polar 
cases of low aerosol and water abundance) , 

3.3.2 Detection of humidity variations via horizontal ranging, 

3.3.3 Upward ranging by means of ground-based Ildars, 

3.3.4 Downward ranging from the 200kM Space Shuttle altitude. 
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These cases are presented not only to discuES the feasibility of Space 
Shuttle observations per se , but also to show the nature of lidar returns - 
and their sensitivity to water vapor concentration - in ejcperiments that 
will naturally be undertaken in the development of a Shuttle lidar capability. 
3.3.1 Uniform Atmospheres (horizontal ranRlng) 

Figures 5 and 6 depict the lidar returns at sea level for two of the 
standard atmospheres adopted by McClatchey et al (1972). Laser energy and 
range resolution are shown on the graphs; detection efficiency is assumed 
to be 10%. As mentioned earlier, the off line returns employ a wavelength 

O 

of 8600 A which is reasonably representative of operation throughout the 

0 

very near infrared (7000 A - ip). The terms "clear" and "hazy" refer to 

standard haze models, whose visibilities are 23 and 5 kM and whose sea 

level aerosol densities are 2.828x10^ and 1.378x10^ cm respectively. 

Water vapor concentrations were calculated at various altitudes by R. G. 

Ellingaon. The "tropical" and "midlatitude summer" values for [H^O] at 

sea level are 6.405x10^^ and 4.719x10^^ cm~^, respectively. These models 

are employed in most of the calculations given here, except for special 

polar or stratospheric cases. Aerosols are usually treated as Mie scatterers 

1 fis 

having a backscatter coefficient 8 -v — f— , Molecules are assumed to be 

TT 5 4v 

Rayleigh scatterers and non-absorbing (except for H 2 O) . 

The figures show both the off line signals and the returns as modified 
by on line operation with different H^O absorption cross sections. The 
ranges of ^^2^^ have been chosen to represent the line strengths avail- 

able in the very near Infrared. The clear, tropical situation offers ad- 
vantages for ground truth experiments, in that the aerosol extinction with 
range is far less drastic, and the higher humidity lends a greater accuracy 
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Figure 5 
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Range (kM) 


Figure 6 
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to the logarithmic derivative which must be calculated from such signals 
in order to extract [H 2 O] (see Equation 6). 

Figures 7 and 8 deal with the opposite extreme of stratospheric 
observations by horizontal ranging from aircraft. Figure 7 compares the 
off line returns to be expected from the molecular atmosphere at 20 IcM 
altitude and a Rayleigh-scattering aerosol whose concentration is 30/cm . 
Variations in an aerosol concentration of this magnitude at this altitude 
will clearly dominate the molecular backscatter contribution. Assuming 
only the molecular off line return, Figure 8 shows the large water vapor 
absorption cross sections required to render 5 ppm of H^O detectable in 
the stratosphere. These H 2 O cross sections can be obtained further into 

O 

the infrared, or at high altitude for the 9400 A lines owing to line- 
narrowing at the greatly reduced ambient pressure in the stratosphere. 

Figures 9 and 10 estimate the lidar returns to be expected at South 
Pole Station (ground level »^3kM altitude) in the summer when there is a 
low altitude, large particle aerosol abundance of order 150 cm (W. Zoller, 
private communication). Mie backscatter is assumed to apply, having a 
magnitude roughly 1/5 that for Rayleigh-scattering particles. Thus the 
aerosol return here is about the same as in the stratospheric case shown in 
Figures 7 and 8, while the molecular backscatter is much larger owing to 
the higher atmospheric pressure. 

The distinct on line return for H20,iu spite of the low value of 
[H 2 O] ««2xl0^^ cm is due to the large absorption cross section assumed, 

O 

which is appropriate for lines in the vicinity of 9400 A (see subsection 3.2), 
We conclude that lines of this strength are mandatory for water vapor assay 
in the polar regions, and that weaker lines at shorter wavelength are not 
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Antarctic Summer (approx.) 
BACKSCATTER COMPONENTS 
horizontal lidar returns 
at South Pole Station 
Aerosol density: ISO cm 
300 mJ, I M rcvr. area 
Iftsec range gate, X~860 nM 
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Figure 10 
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suitable for this type of Ildar operation. 

°g-te,°tlou of Humidity Varletlou e (horlroor.l 
Several Ildar calculation, ere given here having to do „lth the 
"vlelblllty'’ of huuidity variation. In a unlforr. or varying aerosol back- 
ground. Figures 11, 12 and 13 compare three altitude cases In »hlch the 
aerosol distribution 1. unifon. (clear) along tb,. line of sight while the 
humidity Increase. In a well-defined •'bump" between 4 and 9 kH range. 
Figure. 11 and 12. for sea level and 3 kH altitude In a midlatitude summer 
atmosphere, show the effects of a humidity profile which locally achieves 
a tropic al level appropriate to each respective altitude. Figure 13 at 
10 kH assume. 20% and 60% humidity Increases above the background value. 

Aside from the characteristic shift of the off line return, as the 

extinction coefficient diminishes with altitude, the tmln difference In 

the three figures Is the value of the H^O absorption cross section required 

to give a distinct humidity signature In the presence of the variable 

profile. The range of cross sections Is compatible with near Infrared 
operation. 

In the case of Figure 11. the on line signal changes by about 3% 

by the time one la observing the 5 kH range position of the fully developed 

profile. This change 1. about lOx the signal shot noise in the correspondin 
rangfe gate. 

It seem, likely that Ildar operation in the very near Infrared will 
be useful for detecting water vapor profiles connected with developing 

precipitatio n or uPHUln£^fjol^ These will be natural 
objects of study during the development and testing of Ildar systems tor 
Space Shuttle application. 
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Flgure 12 
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Figure 13 





-52- 


Figures 14 - 17 depict the different types of lidar returns one will 
encounter when both the aerosol and water vapor concentrations vary along 
the line of sight. Various combinations of (clear, hazy) and (tropical, 
midlatitude summer) are considered. Figures 14 and 15 show the probing of 
the edge of a large region of atmospheric variation, such as one might 
encounter in horizontal ranging perpendicular to a sea coast. The manifest 
non-parallelism of the on line and off line returns Implies measurability 
of IH 2 O] in. these gradient regions, given sufficient transmitted laser 
energy. Figures 16 and 17 are concerned with mutual and opposite variations 
of aerosols and water vapor, over "bump" portions of the atmosphere which 
are roughly as large as they are distant from the lidar station. The lidar 
signatures are very different, and again' admit of water vapor measurements 
if one employs sufficient laser energy. While the accuracy of these parti- 
cular measurements is not assessed here, it is clear from the variability 
of these signatures between different atmospheric models that many types 
of situations can be distinguished from one another in field operations. 

3.3.3 Vertical Ranging from Ground Stations 
For both this case and the Shuttle case (3.3.4) it was necessary to 
represent the altitude variation of the pertinent optical properties of the 
earth's atmosphere; viz. the extinction and backscatter due to aerosols and 
molecules, and the line absorption due to water vapor. Our emphasis has 
been to establish feasibility rather than to obtain extremely accurate 
simulations of lidar returns. Therefore we have assumed H^O absorption 
cross sections which do not depend on altitude, as well as employing the 
model wavelength of 8600A and the "1/5 X isotropic" assumption for back- 
scatter by Mie-scattering aerosols. The present calculations therefore 



# photons detected /I50M range cell 


HgO LtDAR RETURNS 

horizontal ranging at sea level 
(clear, midlatitude summer} 

— »-(hazy , tropical) 

300 mJ, IM^ rcvr. area 
l/isec range gate, X~860nM 

^ • a C ^ ^Qt )5 ^2*^ 


I2E-24) 


< 

C5.5E-24) 


OFF LINE 


fE-24J 


HgO vapor density 


1 6.405 E 17 cm' 


4.719 E 17 cm 


Aerosol density 


13,780 cm 


2,828 cm 


Range (kM) 


Figure 14 



le cell 


- 54 - 


HgO LIDAR RETURNS 

horizontal ranging at sea level 

{clear, tropical) | 

(hazy, midlat. summer) 

300 mJ , IM^ rcvr. area 
l^sec range gate, X~860 nM 

rTTlT) cr , H,0 (eff.) 10^ cm^ 

V ^ abs ^ 






# photons detected / 150 M range cell 


-55 



,'V\ OF TBE 
VOOE 


Figure 16 


# photons detected / 150 M range cell 


-56- 



Figure 17 


-57- 


err on the pessimistic side as regards the detectability of H 2 O at high 
altitude, where the lowered pressure and temperature in fact enhance the 
absorption cross section at line centers and therefore impart a larger 
optical depth 3^^ to laser radiation tuned to 

Table 5 sh' a the curve - fitting relations which were employed in 
order to carry out the full lidar simulations on a Hewlett-Packard 9810/9862 
calculator/plotter. Information on these programs will be reported in the 
near future and will in the meantime be available on special request. 

Figure 18 shows the total off line return for a ground based lidar 

O 

operating near 1 •= 8600 A, ranging upward into a standard clear, midlatitude 
summer atmosphere. Also, fictitious returns are shown for purely aerosol 
or purely molecular scattering and extinction. Figure 19 then shows the 
corresponding on-H^O-line returns for different absorption cross sections, 
compared to the off line return. Figure 20 shows the ratios (off line/ 
on line) as functions of range and cross section. 

The marked rate of decrease in signal at low altitude is due to the 
preponderance of aerosols there, as contrasted to the slower fall off at 
higher altitude where molecular backscatter is dominant. It is clear that 
upward-ranging lidar experiments will in general present a problem in the 
dynamic range of measurement. Careful optimization of line strength is 
required in order to obtain reasonable accuracy in [H 2 O] above 2 kM; 
sensitivity to the lower humidity there requires operating the lidar on 
lines which show a strong extinction relative to the off line return - 
and which therefore create difficulties in measuring said humidity ac- 


curately. 
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Table 5. Typical Atmospheric Models 
(R ” altitude in kM) 


A, Air Molecular Density Profile 

altitude (kM) 

3 

density (gm/cm ) 

0-9 

, „„„ ,„-0.0459R 

1,220 X 10 

9-25 

1 onn .--0.0662R 

1,900 X 10 

25-30 

13,684 X 

R>30 

682.63 X 


B. Aerosol Density Profile (Clear: 23 kM visibility) 

altitude (kM) 


2-4 

4-6 

6-13 

R>18 


_3 

density (cm ) 


0-2 +65.243R^ - 26.847R^ - 1352. 7R + 2828 cm ^ 


-65.243(R-2)^ + 364.61(R-2)^ - 677.2(R-2) + 537.1 
6.8624(R-4)^ - 26.847(R-4)^ - i.6715(R-4) + 119.2 


76.804 10 


-.013917R 


5256.2 10 


-.11019R 


C. Water Vapor Profile (Midlatitude Summer) 
altitude (kM) 


-3 

density (cm ) 


0-10 

6.2E17 X 10 cm 

10-15 

1.3002E18 X 

15-20 

1.16E14 X 

20-25 

3.0983E12 X io'^‘03445R 
-.0203R 

25-30 

7.3E13 X 10 

R>30 

1.2053E16 X 
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Signal ratio N___/N-,,(H«0), X~700~I000 nM 
OFF ON i. 

Vertical sounding from ground. Clear midlatitude summer 



Figure 20 
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3.3.4 Space Shuttle (doTOward ranging, 200 kM altltu d^ 

By means of the same representations of the atmosphere mentioned in 
3.3,3, we have obtained a number of simulations for a water vapor lidar 
mounted in the Space Shuttle at a (nominal) altitude of 206 kh. Figure 21 
illustrates the basic, off line returns from the lower atmosphere assuming 
the two standard haze models used heretofore. The hazy case shown appears 
to be just about optimum for the maximum lidar return from the bottom of 
the atmosphere. Figure 22 depicts the good sensifvlty to [H^O], using 
appropriate cross sections for the several groups of H 2 O lines in the near 
infrared. For these shuttle cases, it is clear that the problem of photon 
statistics is very pressing and calls for adopting as large a laser energy 
and lidar receiver as possible so as not to sacrifice range resolution or 

accuracy of the humidity measurement. 

Figures 23 and 24 present a broader picture of the Space Shuttle lidar 

for both the hazy and clear midlatitude summer atmospheres. The off line 
return is followed out to an altitude of 30 kM. where only 5 - 6 photons 
will be detected on the average for the nearly 1/2 kM range cell. Above 
10 kM one needs very large cross sections even to detect, let alone measure, 
the water vapor. Therefore we conclude that a nadir-looking lidar measure- 
ment of [H^O] will succeed mainly in the troposhpere - and conceivable very 
well there - and probably not in the stratosphere. The use of many lidar 
pulses is limited by the average laser power permitted in the Shuttle ap- 
plication, and by the desired spatial resolution in the atmosphere. The 
distance scale of interest in the stratosphere is much larger chan at low 
altitude, so it may be feasible to accumulate several lidar returns to 
Improve the photon statistics for observations at 10 kM and above. 


photons detected /450 M range cell 
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Figure 21 
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Flgure 24 
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Figures 25 and 26 for a standard tropical atmosphere complement 
Figures 23 and 24. They demonstrate primarily a slight enhancement of the 
absorption by H 2 O, which would not materially alter the design of such a 
lidar experiment. 

Figure 27 describes the origins of the off line return as seen from 
the Space Shuttle, in terms of Rayleigh scattering from the molecular 
atmosphere and (assigned) Mie scattering from aerosols. If in fact the 
stratospheric aerosols scatter more isotropically than those at low alti- 
tude, a stratospheric signal enhancement over the present numbers will result. 
This, together with increased absorption cross sections owing to line- 
narrowing in the stratosphere, may make the high altitude assay for [H 2 O] 
more tractable than it seems at present. For this reason, the stratospheric 
lidar calculation needs now to be tackled via more elaborate representations 
of the scattering and absorption parameters as functions of altitude. 

Figure 28 represents an approximate treatment of orbital lidar 
soundings over the South Pole, using a simplified atmospheric model to see 
what a Shuttle-based instrument is capable of in cold regions where the 
water vapor density is low. Three off line returns are shown for Shuttle 
altitudes of 100, 150 and 200 kM; the backscatter is due to molecules and 
two aerosol layers whose extents are indicated on the figure. The upward 
Jump near groimd level (3 kM) is not a ground return from snow and ice, 
but arises from low level aerosols below a 600 M inversion layer. Three 
on line returns are shown with the 200 kM curve; the highest of these is 

Q 

appropriate for the 9400 A lines alluded to earlier. 

From these curves it is clear that water vapor assay can be carried 
out from orbit in the polar regions, and that reasonably strong H 2 O lines 
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Figure 28 
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are needed. Since the ground return will be a large signal, total column 
content of H^O constitutes an attractive measureinent under these conditions. 

The importance of knowing [H^O] in the polar atmosphere is due in 
part to the scavenging of airborne, condensible pollutants which probably 
takes place there. The precipitation of water vapor as high altitude ice 
crystals must play a role in any such process. The capability of a 
Shuttle-based lidar to see both the vapor and condensed phases of H^O 
will be important in understanding polar precipitation. There may be no 
other method for wide area studies of this problem in remote regions of 
the earth. 

3.4 Accuracy and Practicality of Orbital Lidar for [H^Q] 

This section summarizes the available information on the accuracy 
with which the concentration of tropospheric water vapor may be measured 
via orbital lidar. Since our simulations do not contain an explicit error 
analysis, we rely on related calculations which are known to yield lidar 
returns comparable in magnitude to those presented here. The main source 
of error is signal shot noise, due to the number of detectable photons 
lying within a given range gate, combined with the logarithmic derivative 
operation (Equation 6) which is required for obtaining material concentra- 
tion from the lidar returns. For a finite rtnge resolution Ar, this is 
equ.ivalent to a double ratio of four signals: the on line and off line 

returns from both the top and bottom of the range cell of interest. 

Remsberg et al (1975) have studied the general performance of dif- 
ferential absorption Ildars, particularly in the favorable wavelength range 

O 

3000 A - 1 micron which includes all the water line groups discussed in 
this report. They conclude: 
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"BeCter than 10% measurements of material concentrations are 
possible below about 15 kM for lOJ/pulse laser energies. In 
the wavelength region 0.3 < A < 1.0 pM, oxygen, water vapor 
and possibly NO 2 can thus be measured in the troposphere." 

2 

This conclusion assumes: Shuttle altitude = 200 kM, telescope area = 1 M , 

overall detection efficiency « 10% (as in the present report), and a height 
resolution of 2 kM. 

V’e can compare this estimate with one scaled from the SRI report by 

Wright et al (1975) on lidar measurements of atmospheric pollutants such 

as NO 2 , using as a scaling factor the available spectral contrast [o(? o(X;] 

-19 2 

(Equation 6) . This is approximately 3 x 10 cm for NO 2 around 4500 A 

—2 2 

according to Wilkerson et al (1974), whereas we may take it as roughly 10 cm 

for H^O lines discussed in the present report which yield distinct on line 

returns (e.g., see Figure 22). Since the SRI estimate of minimum detectable 

[NO 2 ] is equivalent to 0.01 ppm for a single lOJ lidar pulse pair, the re- 

—2 -IQ -22 

lated quantity for H 2 O is 10 ppm x 3 x 10 /lO = 300 ppm. Moreover, 

4 

since the typical H 2 O abundance is of order 10 ppm, the minimum detectable 
variation in [H 2 O] is therefore approximately 300/10^ = 3% under the given 
lidar assumptions. 

Further, Remsberg et al continue: 

"... more detailed simulations were conducted for water vapor 
in the 0.72 pM region using a O.lJ/pulse laser operating at 
1 pps. Relative humidity data accurate to 4% can be obtained 
from a Shuttle platform for a 0.5 kM range cell near the sur- 
face and for horizontal grids of 250 kM." 

These lidar conditions are equivalent to about 4J total energy output during 
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the Shuttle transit over the horizontal grid spacing, or about 1 Joule 
for a 4 X larger vertical range cell of 2 kM. That this laser energy 
estimate for 4Z relative humidity accuracy is about 1/10 of that extrapo- 
lated from the SRI results on [NO^] is encouraging, and may reflect the 

noticeably greater off line scattering extinction at 4500 A as opposed to 
7200 A. 

It is clear that integrated laser pulse energies of order 3J will be 
needed for tropospheric water vapor observations with a 1 telescope in 
the Shuttle, in order to obtain ~1 kM vertical resolutions and relative 
humidity accuracy better than 5%. Since 300 mJ pulses can be expected with 
a 3J pump laser and 10% conversion efficiency, this would require 10 pulses - 
or a 10 receiver in the case of one 300 mJ tunable, near IR pulse. Very 
large receiving telescopes have Indeed been considered for the Shuttle, and 
we would conclude that they should be made available for Ildar measurements 
in order that the maximum useful information on the atmosphere be obtained. 

In conclusion, one can write a criterion for useful near IR Ildar 
measurements (Shuttle altitude - 200 kM) giving relative hLimidities accurate 
to roughly 3% at low altitudes in the atmosphere, namely: 

N*E*A*K*AR 5 3 X 10^ J'M^, ( 

(near IR, N pulses, 3% humidity accuracy) 

where AR is the height resolution (lidar range cell), A is the receiving 
telescope area, K is the overall detection efficiency, E is transmitted 
laser energy per pulse, and N is the number of lidar pulses accumulated 
for a given observation. 
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Another criterion may be written for the average electrical power 
P (watts) required by the lidar apparatus during a time interval of observa- 
tion which covers a horizontal distance D (meters) at ground level: 


n7.D-A.K.AR i 2 X 10^ W-M^, 


( 16 ) 


where n is the conversion efficiency from electrical power into average 
lidar output power . This relationship roughly accounts for the ground 
speed of the Shuttle (^^7 kM/sec) and the estimate (above) that the product 
NE must be of order 3J. 

An example of the consequences of these criteria is important to cite 
in closing this section. For mesoscale observations (Section 2) it is 
preferable to have a horizontal resolution capability of order 1/10 the 
size of 100 kM spatial variations; i.e., set D = 10 kM. Given a detection 
efficiency K ■= 10% and a height resolution AR = 1 kM, this requires 
nP*A “ 2 WM^. For an average electrical power of 1 kW dedicated to the 
Shuttle Lidar Facility, and q “ 10 ^ , the required receiver area is 2 ; 

this is considered to be a reasonable value by most lidar design groups. 

Thus a Shuttle-based water vapor lidar instrument appears to be practi- 
cal. The magnitude of the atmospheric returns and the absorption strengths 
of resolvable H 2 O lines are large enough for useful meteorological studies of 
the altitude distribution of water vapor. Moreover, an orbital lidar probably 
is the only means for worldwide surveys of tropospheric H^O, including many 
remote regions of the earth. 


This figure mediates between the value of 5 x 10 we estimate as rea- 


sonable for proven laser systems and 3 x 10 ^ which may soon be possible 

RK1'K> UF THE 

with excimer-laaer-pumped tunable lasers. 
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4 . Exploratory Laser Experiments 


This section describes experiments undertaken at the University of 
Maryland to explore the application of various dye laser methods for 
generating laser pulses which could be tuned over H 2 O absorption lines in 
the visible and near Infrared. The work reported here is the basis for 
subsequent laser developments under NASA grant NSG - 1156, which have been 
summarized by Mcllrath et al (1975) and Wilkerson et al (1975) and will 
be described in greater detail in the near future. Here we report on the 
following: 

(4.1) Operation of a long pulse dye laser (10 ^ sec full width at 
base) with a tunable, narrow band output at high energies near 

O 

the 5915 A water vapor absorption bands. 

(4.2) Assembly and operation of a short duration dye laser (2 x 10 ^ 

O 

sec full width at base) near the 5915 A water vapor absorption 
bands . 

(4.3) Construction of a dye laser to be pumped by the lasers in 

(1) and (2) above to operate in the red and near infrared. This 

O c 

laser has operated in the region of the 6513 A and the 6981 A 
absorption bands. 

(4.4) Preliminary studies of the beam divergence of the output of 
the laser-pumped systsm. 

In the rest of this section we discuss these results in more detail. They 
are summarized in Section 4.5. 

4.1 Operation of Long Pulse Dye Laser 

The long pulse dye laser is a high energy system capable of 1.25 Joules 
output in a 10 ^ second pulse. We have used it in a cavity with a diffraction 
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grating as the tuning element and have kept the output energies below 
0.4 Joules. No grating damage was observed. Using Rhodamine 6G dissolved 
in methanol we tuned to 5900 A with a line width of < lA. The beam diver- 
gence was studied by photographing the beam with a four lense camera to 
get the intensity distribution at various distances from the output. 
Preliminary measurements indicate that the central core of the beam had a 
divergence of less than 6 x 10 rad. This system is now being used in 
conjunction with a sodium vapor heat pipe oven to tune onto the narrow 
sodium resonance lines at 5890 A, Observations of resonance scattering for 
many shots should give a good indication of shot-to-shot spectral stability 
in this configuration. Tuning through the resonance lines will also give 
a measure of the spectral width. 

^ • 2 Short Pulse Fla shlamp— Pumped Laser 

A flashlamp-pumped dye laser using a Marx Bank capacitor storage 
was put into operation. Outputs of 184 mj were obtained in the 5800 A 
region using Rhodamine 6G in methanol. It would be possible to obtain 250 mj 
of energy by optimizing the dye concentration. The pulse width is less 
than 250 x 10 seconds at the base and 10 seconds at half maximum giving 
an output power of 2 x 10^ watts. This system was also operated with 
cresyl violet to obtain output at longer wavelengths. The output mode 
structure was ■ ompllcated, as expected, and the flashlamp-pumped dye laser 
is used mainly as a pump for the laser-pumped dye laser. Direct output 
of the cresyl violet dye laser was 45 mj with no effort made to optimize 
the output. 
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4.3 Laser-Pumped Dye Laser with Creayl Violet 

A dye laser was constructed consisting of a stainless steel dye cell, 

a total reflector, and an output coupling reflector. The cell is designed 

to be driven by the output of a flashlamp-pumped laser, either dye or ruby. 

The laser was operated with cresyl violet in methanol and pumped by the 

output of the fast pulsed laser operating with Rhodamine 6G. Conversion 

efficiencies of 4% were achieved without efforts at optimizing. The output 

of cresyl violet is near 6500 A. Photographs of the output beam show a 

-3 

beam divergence of less than 2 x 10 rad. 

4.4 Laser-Pumped Dye Laser with Nile Blue Perchlorate 

The laser-pumped dye laser was operated with a solution of Nile Blue 
Perchlorate in Methanol. This system lased well with an output of 1.2 mj 

O 

centered on 6932 A. The output should be increased by changing the pump 
wavelength and dye solvent and by optimizing the dye concentration. The 
mode pattern is similar to the cresyl violet case. 

4.5 Stinnnary 

Enough information has been gathered to allow construction of a tunable 
dye laser operating in the near infrared region which is suitable for 
differential absorption Lidar experiments. Water vapor absorption lines 

O O 

in the 7150 - 7400 A region and the 8100 - 8400 A region offer a wide 
variety of line strengths and can easily be reached with tunable dye lasers 

Q 

(see Section 3.2). The H2O absorption lines are approximately .07 A in 
width, a width which again is easily matched by dye lasers. The next step 
in this work has been the construction of such a tunable dye laser that is 
adequate for prototype experiments. 
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For ease of construction, a dye laser pumped by a giant pulsed ruby 
laser is most attractive. Experiments on such a laser are currently under- 
way under NSG-1156. Sufficiently energetic and spectrally pure pulses are 
now available for "ground truth" measurements of atmospheric water vapor, 
according to simulations of the kind represented by Figure 19 . The ease 
of reaching this experimental configuration assures us that the laser 
parameters needed for Shuttle observations of H^O can readily be obtained 
with space-qualified instruments during the late 1970' s. This further 

buttresses the case for the practicality of orbital H2O lidar from the 
Space Shuttle platform. 
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Appendlx 1 

The relationship between different meteorological 
parameters used to specify the water 
vapor content of air 

It is useful to briefly discuss the paramatere used by meteoro- 
logists to denote the water vapor content of air. The most commonly 
used parameters are: 

1. Mixing ratio oj H P^/P^ 
p^- water vapor density 

p,- Density of dry air containing the water vapor 
d 

2. Specific humidity q = P^/Pg 

p - density of moist air containing the water vapor 
q ” w/(l + m) 

3. Percent Relative humidity r = 100 e/e^(T) 
e - water vapor partial pressure 

e - saturation vapor pressure with respect to a plane liquid 
water surface at the temperature, T, of the vapor. 

4. Absolute humidity p^ 

It should be noted that ui and q are conservative properties of air 
undergoing dry adiabatic processes. Thus, these parameters serve as 
convenient tracers of air for many atmospheric problems. 

The absolute humidity or water vapor density is used less than m, 
q or r. However, or the number density is the variable obtained 
from the Ildar retrieval technique. The water vapor data requirements 

; '-M' STY OF THE 
ORKJIYAL PAOF IS POOR 
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dlscusaed refer to relative humidity. The material below relates relative 
humidity to water vapor density. 

i'3suming water vapor obeys the ideal gas law, one may relate 
r to by the following equation. 


r 


100 p 


R T 
w 

V e (T) 
8 


where is the gas constant for water vapor. 

^®t Tg, p^^ and Tq be the correct values of r, p^ and T, respec- 
tively, and let 6r, 6p^ and 6T be the uncertainty of r, p^ and T, res- 
pectively. Assuming the Clausius-Clapeyson equation is approximately 
correct for finite temperature differences, then the fractional error 
of relative humidity can be related to the fractional error in density 
and temperature by the equation 


r 

o 


<5P. 


+ (1 - 


6T 


VO 




where is the latent heat of vaporization. 

It is easily seen that the relative error in r equals the relative 
error in p^ when the temperature is known precisely. However, if the 
temperature uncertainty gets large, it will mask the effects of the water 
vapor uncertainty. 


- 0.688 


L 

At T 0“C, Che term (1 - 

^ o 

for 6T - 1.0, 5.0 and 10.0'’C, 


■) amounts 
^o 

respectively. 


to -0.063, -0.344, 
This effect of the 


temperature uncertainly should not be overlooked when estimating the un- 
certainty in relative humidities obtained from water vapor density estimates. 
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